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Abstract 

We provide the representation of quasi-periodic solutions of periodic Dirac 
equations in terms of the spectral parameter power series (SPPS) recently 
introduced by V.V. Kravchenko [lj[2]|3]. We also give the SPPS form of 
the Dirac Hill discriminant under the Darboux nodeless transformation 
using the SPPS form of the discriminant and apply the results to one of 
Razavy's quasi-exactly solvable periodic potentials. 

1 Introduction 

The connections between the Dirac equation and the Schrodinger equation are 
known since a long time ago ^4 and have been strengthen in the supersynimetric 
context soon after the advent of supersymmetric quantum mechanics in 1981 
[51 El [Zl IS]- There are currently interesting applications of this approach in 
condensed matter physics [HI dHl El- In this work, we are interested in the 
same connection in the case of periodic potentials, see e.g. [H]. We here write 
the Dirac Bloch solutions in Kravchenko form (power series in the spectral 
parameter) and also the Dirac Hill discriminant in the same form and apply the 
results to an interesting quasi-exactly solvable periodic potential. 



* Operator Theory: Advances and Applications 220, 153-162 (2012). Springer, Basel AG. 
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2 Schrodinger equations of Hill type 



The Schrodinger differential equation 



L [/(x, A)] 



(1) 



with T-periodic real- valued potential q[x) assumed herewith a continuous bounded 
function and A a real parameter is known as of Hill type. We begin by recalling 
some necessary definitions and basic properties associated with the equation ((T)) 
from the Floquet (Bloch) theory. For more details see, e.g., [T3l [14] . 

For each A there exists a fundamental system of solutions, i.e., two linearly 
independent solutions of ([l}, A) and f2{x, A), which satisfy the initial con- 
ditions 



Then the Hill discriminant associated with equation ([ij is defined as a function 
of A as follows 



The importance of D{X) stems from the easiness of describing the spectrum of 
the corresponding equation by its means, namely |13j : 

(1) sets {Ai} for which |-D(A)| < 2 form the allowed bands or stability 
intervals, 

(2) sets {Aj} for which |-D(A)| > 2 form the forbidden bands or 
instability intervals, 

(3) sets {Afc} for which |I?(A)| — 2 form the band edges and represent 
the discrete part of the spectrum. 

Furthermore, when D{X) = 2 equation ^ has a periodic solution with the 
period T and when D{X) = —2 it has an aperiodic solution, i.e. f{x -I- T) = 
— /(x). The eigenvalues A„, n = 0, 1, 2, ... form an infinite sequence Aq < Ai < 
A2 < A3 < A4..., and an important property of the minimal eigenvalue Aq is 
the existence of a corresponding periodic nodeless solution u{x,\q) [13]. The 
solutions of (H]) are not periodic in general, and one of the important tasks is 
the construction of quasiperiodic solutions defined by f±{x + T) — /3±(A)/±(x). 
Here, we use James' matching procedure |15j that employs the fundamental 
system of solutions, /i(a;. A) and f2{x, A), in the construction of the quasiperiodic 
solutions as follows 



/i(0,A) = l, /((0,A)=0, /2(0,A)=0, /^(0,A) = 1. 



(2) 



i^(A) =/i(T,A) + /^(T,A). 



/±(x. A) = (3l{\) [h{x ~ nT, A) + a±/2(x - nT, A)] 




(3) 



where a± are given by |15) 



a± = 



/^(T, \)- h{T,\)T{D^ 
2/2 (T, A) 



(A)~4)^ 



(4) 
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The Bloch factors (3±{X) are a measure of the rate of increase (or decrease) in 
magnitude of the Unear combination of the fundamental system when one goes 
from the left end of the cell to the right end, i.e., 

. /i(r,A)+a±/2(r,A) 
' /i(0,A)+a±/2(0,A) ■ 

The values of (3±{\) are directly related to the Hill discriminant, /3± (A) = 
^{D{X) =F y/D^{X) - 4), and obviously at the band edges ^+ = /3_ = ±1 for 
D{X) = ±2, respectively. 



3 SPPS representation for solutions of the one- 
dimensional Dirac equation 

We consider the following Dirac equation 

L[W]^[-iayd^ + (T^<^{x)]W = ujW, (5) 

where the scalar potential ^(x) is periodic function with period T, W is the 

spinor ^ ^ ^ ^ ^ ^^''^ '^^ Pauli matrices '^^ ~ q ^ ^'^"^ 

-i 

1 

The uncoupled Schrodinger equations derived from equation ([5]) are 

M, + $)(d, + $)/ = A/, (6) 
(d, + <i>)(-d, + <i>)g = A5, (7) 

where A = is the spectral parameter. It is clear that the solutions / and g 
are related by the following relationship 

(d, + $)/ = wg, (8) 

therefore with the solution / at hand, we can construct the solution g immedi- 
ately. 

We start with equation ([6]). Notice that the solution u of the equation ([6]) 
for A = can be obtained as follows u{x) = ^ '^^^^'^^ and u{x) is a nodeless 
periodic function with the period T if ^{x) e and ^{x)dx = 0. 

Once having the function u{x) the solutions fi{x, A) and f2{x, A) of ([6]), ([2]) 
for all values of the parameter A can be given using the SPPS method [1] . 

h{x,X) = 4^Eo(x,A)+«'(0)«(a;)I]i(x,A), 

(9) 

f2ix,X) = — u(0)u(a;)Si(x, A). 
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The functions Eq and Si are the spectral parameter power series 

oc oo 

So(x, A) = ^ x(2")(a;)A", A) = ^ X^^"-^\x)X"-^ , 

n— n— 1 

where the coefficients X(")(a;), are given by the following recursive 

relations 

= 1, = 1, 

(JoX^"-^'>{CW{^)d^ foranoddn 
= i (10) 
[-/;i("-i)(0^ foranevenn 

r-/„"x("-i)(C)^ foranoddn 
X(-\x)^l (11) 
[/p'x("-i)(0w2(^)^^ foranevenn. 

One can check by a straightforward calculation that the solutions /i and /2 
fulfill the initial conditions ([U, for this the following relations are useful 

(So(x,A))' =-^4^, whereEi(x,A) = X(^^^-'\x)X^ (12) 
u x — : 



n=l 



and 



(Ei(x,A))l = -^4rr' ^l^ere Eo(a;,A) - X(2")(x)A". (13) 

^ ^"^^ n=0 

The pair of linearly independent solutions gi{x,X) and (72(2;, A) of ([7]) can be 
obtained directly from the solutions (|9]) by means of ([8]). We additionally take 
the linear combinations in order that the solutions gi{x, A) and 172(2^1 A) satisfy 
the initial conditions gi{0, A) = (72(0' A) = 1 and 5^(0, A) = .92(0, A) = 

51 (x, A) = ^Eo(x, A) - f f M x,\), 

U[X) AU(U)U(X) 



(14) 



Xu[0)u[x) 

Thus, the two spinor solutions of the Dirac equation ([5]) are given by 

Wi = ( •^M , and W^2 = I 
and these solutions satisfy the following initial conditions 
W^i(O)- (I), andM^2(0)= ( J 
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3.1 Bloch solutions and Hill's discriminant 



The second order differential equations ^ and © have periodic potentials 
Vi,2 = <i>^ =F ^' , correspondingly. The important tasks for this case are the 
construction of the Bloch solutions which are subject to the Bloch condition 
f{x + T) = e**''"^/(a;) (fc, a wave number) and the description of the spectrum. 

In [TB] the SPPS representations of Hill discriminants £'/(A) and Dg{X) 
associated with the equations © and ([7]) were obtained in the form 

I?,(A) = ^So(T,A) + 4?^So(T,A) 

+ (u'(O)u(T) - u{0)u'{T)) Si(r, A) , 



n^W - ^^o(T, A) + ^£o(T, A) + iu'iOHT)- 
-u'(r)«(0))Si(T,A). 

It is clear that since u{x) is a T-periodic function (u(0) = u{T) ) the expression 
in brackete in the above formulae vanishes. Now writing the explicit expres- 
sions for Eo(T', A, Ao) and Eo(r, A, Aq), a representation for Hill's discriminant 
associated with ([6]) and ([7]) is the following 

oo 

Df{X) = Dg{X) = ^ (^X^^^^T) + x(2")(t)) A". (15) 

ri=0 

Equations ^ and Q are isospectral and we obtain the Hill discriminant asso- 
ciated with the Dirac equation ([5]). We formulate this result as the following 
theorem: 



Theorem. Let $(x) G be a T -periodic function which satisfies the condition 
Jp ^{x)dx = 0. Then the Hill discriminant for ^ has the form 

oo 

DH'(w)=5](^^'"'(r)+X(2«)(r)) W^n ^ 
n=0 

where X^'^"-') and X^'^"''^ are calculated according to i fT^) and i f77]j . u — e"^*^^)'^^ 
and the series converges uniformly on any compact set of values ofuj. 

In order to construct the Bloch solutions for the Dirac equation ([5]) we use 
the solutions (|9]) and ([T4l) and apply the procedure of James [TF. Notice that 
because the Hill discriminants for the equations © and ^ are identical the 
Bloch factors for both equations are equal. The so-called self-matching solutions 
for the equations (O and ([7]), are correspondingly 

F±{x,\) = A) +a±/2(a:;, A) and G±(x, A) = gi{x,\) + ^±.92(2;, A), 
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where a± and b± are calculated by the formula (|31) with the corresponding 
fundamental system of solutions ^ and (IT4l) . By means of F± and G± we 
write the self-matching spinor solution of the equation ([5]) 

Finally, the Bloch solutions of the equation ([5]) take the form 

W±ix,\)=(3l{X){w±{x-nT,\)), I "Jfo!±l,±2t!.^^ 

4 Numerical calculation of eigenvalues based on 
the SPPS form of Hill's discriminant 

As is well known [T3], the zeros of the functions D{X) =F 2 represent eigenval- 
ues of the corresponding Hill operator with periodic and aperiodic boundary 
conditions, respectively. In this section, we show that besides other possible 
applications the representation (jlSp gives us an efficient tool for the calculation 
of the discrete spectrum of a periodic Dirac operator. 

The first step of the numerical realization of the method consists in calcula- 
tion of the functions X*^"^ and X^") given by (ITUl) and (ITT|) . respectively. This 
construction is based on the eigenfunction u{x). Next, by truncating the infinite 
series for D{X) (fTSl) we obtain a polynomial in A 

N 



Dn{X) = + X(2")(r)^ A" (16) 

n=0 

N 

= 2 + ^ + x(2»)(r)) A". 



The roots of the polynomials I? at (A) =F 2 give us the eigenvalues corresponding 
to equation ([1]) with periodic and aperiodic boundary conditions, respectively. 

As an example we consider the Dirac equation ([S]) with the scalar potential 

^ ' [2 i~A{i)co%2x 

with j4(^) = ^1 — + and ^ a real positive parameter. This scalar po- 
tential satisfies the conditions of theorem 3.1. The corresponding second order 
differential equations are 

-£-^g + V2g = A.g, 
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where the Schrodinger potential 



14(a;) = ^(l-cos4a;)-3Ccos2a; , (17) 
8 

is the case m — 2 in the quasi-exactly solvable family of the so-called trigono- 
metric Razavy potentials [17], Vr = ^ (1 — cos 4a;) — (m + 1)^ cos 2a;. For a 
given integer m, if ^ < 2(rn + 1) the potentials Vj^{x) are of single-well periodic 
type and if ^ > 2(to + 1) they are of double-well periodic type. 



y.(.) = M^)+4cos2.(i- (18) 



is the supersymmetric partner potential and therefore it is also quasi-exactly 
solvable. The Schrodinger equations with these potentials can be used for the 
description of torsional oscillations of certain molecules [17 . Plots of the po- 
tentials Vi{x) and V2{x) are displayed in Fig. 1 for two values of 




Fig. 1: The Razavy potentials Vi (solid lines) for £, — 2 and ^ = 11 given by fl7]) and 
its partner potentials V2 (dashed lines) as given by for the same values of ^. 



The computer algorithm was implemented in Matlab 2006. The recursive 
integration required for the construction of Xq^\ Xo^\ X^"^ and X^"-^ was done 
by representing the integrand through a cubic spline using the spapi routine with 
a division of the interval [0,7r] into 5000 subintervals and integrating using the 
fnint routine. Next, the zeros of Z3jv(A) ± 2 were calculated by means of the 
fnzeros routine. 

In the following tables, the eigenvalues were calculated employing the SPPS 
representation for four different values of the parameter ^. The first two 
values are below the threshold value ^thr = 6 for m = 2 from single-well to 
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double-well types of Razavy's potentials while the last two values are above this 
threshold value. For comparison, we use the eigenvalues given analytically by 
Razavy in terms of the parameter ^ as follows [T7| 

Ao = 2 (l - yi+e^) , A3 = 4, A4 = 2(1 + ^1+^2) • 





e = i 






n 


A„ (SPPS ) 


A„ (Ref. [T7I ) 







-0.828427124746190 


-0.828427124746190 




1 


-0.628906956748252 






2 


2.315132548422588 






3 


3.999991462865745 


4 




4 


4.828420096225068 


4.828427124746190 




5 


9.238264469324272 






6 


9.294265517212145 








e = 2 


e = 2 




n 


A„ (SPPS ) 


A„ (Ref. El) 







-2.472135954999580 


-2.472135954999580 




1 


-2.428136886851045 






2 


3.184130151531468 






3 


4.000004180961838 


4 




4 


6.472138385406806 


6.472135954999580 




5 


9.864117523158974 






6 


10.253256926576858 








e = ll 


? = 11 


n 


A„ (SPPS ) 


A„ (Ref. dl]) 





-20.090722034374522 


-20.090722034374522 


1 


-20.090721031408926 




2 


3.999728397824670 




3 


4.000000543012631 


4 


4 


24.092379855485746 


24.090722034374522 


5 


24.125593160436161 




6 


36.212102534969766 






C = 20 


^ = 20 


n 


A„ (SPPS) 


A„ (Ref. [T71) 





-38.049968789001575 


-38.049968789001575 


1 


-38.049968788934475 




2 


3.999999942823312 




3 


3.999999999630503 


4 


4 


42.050313148383374 


42.049968789001575 


5 


42.050347742353317 




6 


74.691604620863302 





In Fig. 2, we display the plots of the Hill discriminants for the values of the 
Razavy parameter ^ = 1, ^ = 2, and ^ = 3, respectively. In general, these 
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plots contain damped oscillations with higher amplitudes at higher ^. On the 
other hand, getting the spectrum in A is equivalent with having the eigenvalues 
uJn — ±VA^ of the Dirac system under consideration. 




-4 -2 2 4 6 8 10 12 14 

A 



Fig. 2: The polynomial Z)jv(A) for the Hill equations with Razavy's partner potentials 
for three values of the parameter ^ calculated by means of formula i|J6'|) for N = 100. 
The first minimum of the discriminant in this plot, i.e., for ^ = 3, goes down to -260.9 
at \ = -2.469. 



5 Conclusions 

In summary, in this work we presented the SPPS form of the quasi-periodic 
(Bloch) solutions of periodic one-dimensional Dirac operators as well as of the 
Hill discriminant. We applied the obtained results to the Dirac system with the 
periodic scalar potential that leads to one of Razavy's quasi-exactly solvable 
periodic potentials. 
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